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Abstract 

Increasing the pressure and temperature 
of the engines of new generation of supersonic 
airliners increases the emissions of nitrogen 
oxides (NO x ) to a level that would have an 
adverse impact on the Earth’s protective 
ozone layer. In the process of evolving the 
implementing low emissions combustor tech- 
nologies, NASA Lewis Research Center has 
pursued a combustion analysis program to 
guide combustor design processes, to identify 
potential concepts of greatest promise, and to 
optimize them at low cost, with short turn- 
around time. The approach is to upgrade and 
apply advanced computer programs for gas 
turbine applications. Efforts have been made 
in further improving the code capabilities of 
modeling the physics and the numerical 
method of solution. Then test cases and 
measurements from experiments are used for 
code validation. 

To provide insight into the combustion 
process and combustor design, two- 
dimensional and three-dimensional codes such 
as KIVA-II and LeRC3D have been used. 
These codes are operational and calculations 
have been performed to guide low emissions 
combustion experiments. 

Introduction 

No x emissions from current technology 
combustors are produced in the near stoichio- 
metric primary zone at a rate that increases 
with flame temperature and residence time. 

As new engine cycles increase combustor inlet 
pressure and temperature, the NO x emissions 
increase correspondingly. Most low NO x com- 
bustor design concepts reduce NOx emissions 
by reducing flame temperature and oxygen 


availability for NOx formation. This can be 
accomplished by burning either fuel lean or 
fuel rich in the primary zone, thus avoiding 
high NOx formation at near stoichiometric 
conditions. 

In the process of evolving and imple- 
menting low emission combustor technologies, 
NASA Lewis has pursued a combustion analy- 
sis code program to guide combustor design 
processes, to identify potential concepts of 
greatest promise, and to optimize them at low 
cost, with short turn-around time. The com- 
putational analyses are evaluated at actual 
engine operating conditions. The approach is 
to upgrade and apply advanced computer pro- 
grams for gas turbine applications. Efforts 
have been made in further improving the code 
capabilities for modeling in the physics and 
the numerical method of solution. The test 
cases and measurements from experiments are 
used for code validation. 

Due to schedule constraints, the High 
Speed Research (HSR) analytical research pro- 
gram is being conducted over a period of a 
little over 5 years as shown in Fig. 1 and 
involves the major milestones. The first mile- 
stone was accomplished with the development 
and use of two-dimensional and three- 
dimensional codes, KIVA-II and LeRC3D, to 
guide low emissions combustion experiments. 
These codes will be updated based on results 
obtained from combustion concept experi- 
ments by the end of fiscal year 1995. 

The overall combustion analytical codes 
evolution plan involves in-house research and 
contracts and grants; and it provides strong 
collaborative relationships and technology 
transfer between industry, universities, and 
government agencies. Figure 2 lists the 


activities and the key participants for the 
HSR Combustion Analytical Research 
Program. 

The remainder of the report describes in 
greater detail a description of computer codes 
used and representative results on low emis- 
sion combustor calculations. 

Description of Computer Codes 

To provide insight into the combustion 
process and combustor design, KIVA-II and 
LeRC3D have been used. These codes are 
operational and calculations have been per- 
formed to guide low emissions combustion 
experiments. 

Description of Computer Code KIVA-II 

KIVA-II (Ref. 1) developed by Los 
Alamos National Laboratory, is one of the 
most developed and validated codes of the 
available for multi-dimensional prediction of 
the in-cylinder combustion dynamics in inter- 
nal combustion engines. There are features of 
KIVA-II that make it well suited for other 
applications, so KIVA-II has been adapted for 
gas turbine combustor applications. 

In terms of modeling the physics, the 
major features of KIVA-II are as follows (see 
Fig. 3): 

• KIVA-II is two- and three-dimensional 
turbulent compressible flow solver of reacting 
multi-component gas mixture with liquid 
spray using an Eulerian-Lagrangian approach. 

• Turbulence is modeled using the k-e 
model. 

• Combustion is modeled by a chemical- 
kinetics-controlled model using global or 
detailed chemical reactions (Ref. 2) or by a 
mixing-controlled model (Ref. 3). The user 
can conveniently provide a chemical kinetics 
mechanism by making appropriate modifica- 
tion to the input data file (Ref. 2). 


• The extended Zeldovich NO x mecha- 
nism is included. 

• Stochastic particle spray model 
includes vaporization, coalescence, and 
breakup. 

• A soot formation/oxidation and a 
radiative heat transfer model are also 
included: 

In terms of numerics, KIVA-II is based 
on the following (see Fig. 3): 

• Time-dependent finite-difference code 
with arbitrary mesh capability, using an 
implicit continuous Eulerian technique with 
conjugate residual iteration for the flow 
solver. 

Description of Computer Code LeRC3D 

LeRC3D is highly advanced code for 
gas turbine combustor applications. LeRC3D 
was developed by Carnegie Mellon University 
with the collaboration/sponsorship of Lewis. 

In terms of modelling the physics, the major 
capabilities of LeRC3D (Ref. 4) are as follows 
(see Fig. 4): 

• LeRC3D is two- and three- 
dimensional code that solves the N-S equa- 
tions for turbulent compressible flow of react- 
ing multi-component gas mixture with liquid 
spray using an Eulerian-Lagrangian approach. 

• The turbulence is modeled by using a 
k-e turbulence model with wall functions, or 
by using low Reynolds number k-e model of 
Chen and Patel, or by using a RNG-based 
k-e model. 

• Modeling of combustion is done by 
two different models: the chemical-kinetics- 
controlled model governed by global or 
detailed chemical kinetics mechanisms of 
hydrocarbon combustion, and the mixing- 
controlled model of Magnussen and Hjertager. 
The user can conveniently provide a chemical 
kinetics mechanism by making appropriate 
modification to the input data file. 
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• The chemical kinetics model used to 
study NO x is the Zeldovich mechanism. 

• The spray model includes the fuel 
vaporization model of Raju and Sirinano. 

In terms of the numerical method of 
solution, LeRC3D is based on the following 
(see Fig. 4): 

• The flow algorithm is a finite-volume 
LU algorithm utilizing van Leer flux-vector 
splitting with the HOPE algorithm of Liou 
and Steffan. Source terms are treated implic- 
itly using Shih and Chyu method, diffusion 
terms are treated implicitly using the proce- 
dure of Shih and Steinthorsson. 

• A grid system is generated by using 
an algebraic grid generation method based on 
transfinite interpolation. 

Results and Discussion 
Fuel Injection and Mixing 

Figure 5 shows the prediction of a swirl- 
ing fuel spray in a nonreacting airstream using 
the KIVA-II code and a comparison with the 
experimental results (Ref. 5). The model air- 
assist atomizer embodies a nonswirl inner 
airstream and a swirling outer airstream 
which help to atomize and distribute fuel 
injected from a core tube. Predicted and 
measured air azimuthal velocity profiles are in 
good agreement. Fuel injection and mixing 
become increasingly important as more air is 
used for the combustion process. Detailed 
models of the interaction of the swirling air 
and the fuel spray can provide valuable 
insight into the effect of different variables 
that presently can only be evaluated experi- 
mentally on a global scale at laboratory test 
conditions. 

As an integral part of the fuel nozzle 
calculations, two-dimensional and three- 
dimensional analysis (Refs. 6 and 7) were 
done to provide the swirling air profiles 
through the swirlers of the airblast fuel nozzle. 
The swirl vane cascade analysis provides inlet 


air profiles to the combustor analyses. A 
systematic design procedure of advanced air- 
blast fuel nozzle swirl vane assemblies resulted 
in low-loss, high-turning swirl vane designs 
were done. Figure 6(a) shows the velocity 
profile just downstream of the nozzle inlet. 
Note the velocity is high with a large swirling 
component. This is required for good fuel 
atomization. Figure 6(b) shows the velocity 
field further downstream in the combustor. 
The recirculation zone is evident by the nega- 
tive axial velocity component. 

Low NO x Combustion Analyses 

Figure 7 summarizes the comparison of 
experimental NOx and CO emission index of a 
lean premixed prevaporized (LPP) burner 
(Ref. 9) with KIVA-II code predictions. The 
simplified kinetics mechanism (Ref. 2) was 
used. The predictions agree very well with 
the test data over the range of equivalence 
ratio and residence time reported. Calcula- 
tions using KIVA-II have been performed to 
guide current low NOx combustion 
experiments. 

KIVA-II was used to perform two- 
dimensional analysis of the Rich Burn/Quick 
Quench/Lean Burn (RQL) combustor to 
provide detailed information on the combustor 
internal flow fields, fuel-air mixing, combustor 
emissions, gas temperature distribution, and 
pattern factor. A two-dimensional axisym- 
metric model was used with propane and pri- 
mary air injected at the inlet of the rich burn 
section and quick quench air supplied to the 
two-dimensional slot in the quick quench sec- 
tion. The upper half above the combustor 
centerline showing the gas temperature profile 
is shown in Fig. 8. The rich zone and lean 
zone equivalence ratios were set at 1.4 and 
0.65, respectively. The gas temperature con- 
tours show core hot gas regions that occur in 
the rich zone and lean flame fronts. This 
figure also indicates that the penetration of 
the two-dimensional jet reaches near the com- 
bustor centerline and that thermal quenching 
occurs in the quick quench section. 
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Three-dimensional analyses of the RQL 
combustor were performed using KIVA-II and 
the LeRC3D. Figure 9 shows velocity vectors 
typical of the quick quench-tenure section 
examined at an axial plane at a downstream 
location of half a slot length from the slot 
centerline (Ref. 10). There are twelve 45° 
slanted slots located symmetrically around the 
perimeters of the quick quench section. The 
slots were rectangular in shape with a radius 
at both ends and with an aspect ratio of four 
and the momentum flux ration considered for 
this case is 60. The quick quench section was 
sandwiched between the converging section of 
the rich burner exit and the diverging section 
of the lean burner inlet. One of the most 
striking physical features of quick quench jet 
injected normally into the mainstream from 
45° slanted slots in the appearance of the 
swirling flow field shown in Fig. 9. This swirl- 
ing flow field is responsible for the quick mix- 
ing required for low NO x formation in the 
quick quench section of the RQL combustor 
calculated using KIVA-II. Figure 10 shows 
the velocity vectors and a strong central recir- 
culation zone downstream of the airblast fuel 
nozzle. Calculations using LeRC3D have 
been performed to guide low NO x combustion 
experiments. 

Concluding Remarks 

As physical submodels and numerics 
accuracy improve, analytical models play a 
vital role in the design and development proc- 
ess of low emissions combustors. In spite of 
deficiencies such as accurate chemical kinetics 
modeling of hydrocarbon combustion, model- 
ing of turbulence and turbulence combustion 
interaction, etc. in the current models, the 
computations that were performed using 
KIVA-II and LeRC3D show that these codes 
are capable of producing reasonable results. 
The qualitative insight and trends result from 
parametric variation of combustor design vari- 
ables provided by these codes prediction is 
considered to be a prime factor for reducing 
development iterations and guiding low emis- 
sion combustion experiments. The lack of 


more detailed experimental data has precluded 
more detailed validation of the codes. As 
reliable experimental data becomes available, 
further computations will reveal the capabili- 
ties and limitations of these codes and the 
modifications necessary to overcome the 
limitations. 
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Combustion analysis and prediction codes 



Initial 20 and 3D 

analyses codes 


Code update 

based on 

experiments 


Prediction 
verification 
with test data 


Combustion concept experiments 


NO x destruction Flame tube NO x Mixing and vaporization; 
additives formation/control lean and rich combustion 


f Low-emission component configuration definition and development) 

i “Tf X 

Concept Concepts selected Combustor rig demo's 

assessment/ for combustor (15-20 El) (3-8 El) 

screening development 

Figure 1. — HSR combustion analytical research key milestones. 



Figure 2. — HSR combustion analytical research organization and activities listing. 
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Physics 

Turbulent compressible flow of 
reacting multicomponent gas 
mixture with liquid spray 

k-e turbulence model with wall 
functions 

Combustion models: 

Chemical kinetics controlled, 
mixing controlled model 

NO x formation model: 

Extended Zeldovich mechanism 

Spray model: 

Stochastic model, vaporization, 
coalescence, breakup 

Soot formation/oxidation 

Radiation heat transfer 


Numerical method 

2D or 3D time-dependent finite 
difference code 

Arbitrary mesh 

ICE method with conjugate 
residual iteration 

Optimal quasi-second-order 
upwind convection 


Figure 3. — Key features of computer code KIVA-EL 


Physics 

Turbulent compressible flow of 
reacting multicomponent gas 
mixture with liquid spray 

Turbulence models: 

K-e turbulence model with wall 
functions 

Low Reynolds no. K-e model 
(Chen & Patel) 

RNG-based K-e model 

Combustion models: 

Chemical kinetics controlled 
Mixing controlled 

NO x formation model: 

Zeldovich mechanism 

Spray model: 

Lagrangian model of Raju & 
Sirignano 


Numerical method 

Grid generation: 

Algebraic method using 
transfinite interpolation 

Flow algorithm: 
Finite-volume, LU, implicit 

Code: 

Efficient and robust 


Figure 4. — Key features of computer code LeRC3D. 
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KIVA azimuthal velocity comparison 
Calc = 070889, vte = 180,30 wte = 22,22; k-e 




Radius, cm 


Experimental rig 


Figure 5. — Comparison of predicted and measured air azimuthal velocity of an experi- 
mental air assist atomizer. 


KIVA-P/H V.Profiles no com.(0, 10, 20%SW) 
Axial location, z = 3.90 cm; 750 cycles 



CO 



(b) KIVA-II downstream. 
Figure 6.— Velocity profiles. 
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N0 X EMISSION INDEX (g NOj /kg CjH .) 



L IGNITER 




Figure 7.—Comparison of experimental and predicted NO x and CO emission index of a lean premixed prevaporized burner. 
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til. 


Air 


. Fuel 
(gaseous 
propane) 


Rich Lean 

burn Quick burn 

zone quench zone 

(0=1.4) zone (0=0.65) 


Simplified combustor model 



Predicted temperature contours 

Figure 8. — Predicted gas temperature profile of a rich burn/quick quench/lean burn combustor. 
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Figure 9. — Predicted velocity vector profile at an axial plane at a downstream location of half a slot length from the slot 
centerline of a quick quench design of the rich burn/quick quench/lean burn combustor. 



Figure 10. — Predicted rich burner flow field characteristics of a rich burn/quick quench/lean burn combustor. 
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